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The properties of galactic cosmic rays are investigated with the KASCADE-Grande experiment in the energy range between 10 
an d 10 18 eV. Recent results are discussed. They concern mainly the all-particle energy spectrum and the elemental composition of 
cosmic rays. 
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1. Introduction 

To reveal the origin of galactic cosmic rays is the main 
objective of the KASCADE and KASCADE-Grande exper- 
iments. The results of KASCADE contributed significantly 
to the understanding of the origin of the knee in the all- 
particle energy spectrum of cosmic rays at energies around 
4 • 10 15 eV. It could be shown that the knee is caused by a 
fall-off in the flux of light nuclei |H2j . Energy spectra for five 
dominant elemental groups could be reconstructed (p, He, 
CNO, Si, and Fe), they exhibit a fall-off behavior approxi- 



mately proportional to the nuclear charge of the elemental 
groups. In the energy region between 10 17 and 10 18 eV a 
transition from a galactic to an extra-galactic origin of cos- 
mic rays is expected [3]. Main focus of KASCADE-Grande 
is to understand the end of the galactic cosmic-ray com- 
ponent through detailed investigations of the energy spec- 
trum and the mass composition in this energy region. 

KASCADE-Grande comprises 37 detector stations 
equipped with plastic scintillators, covering an area of 
0.5 km 2 to measure the electromagnetic shower component 
[I]. It also includes the original KASCADE experiment, 
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Fig. 1. Measured lateral distribution of a single event for charged 
particles and muons [9]. 

consisting of several detector systems [5] . A 200 x 200 m 2 
array of 252 detector stations, equipped with scintilla- 
tion counters, measures the electromagnetic and, below 
a lead/iron shielding, the muonic parts of air showers. 
A 130 m 2 streamer tube detector, shielded by a soil-iron 
absorber serves to reconstruct the tracks of high-energy 
(E^ > 0.8 GeV) muons 6J. An iron sampling calorimeter 
of 16 x 20 m 2 area detects hadronic particles [7J. It has 
been calibrated with a test beam at the SPS at CERN up 
to 350 GeV particle energy [8] . 

As an example for a measured air shower, the lateral 
distributions for charged particles and for muons are shown 
in Fig.[T][9j. The charged-particle density is measured with 
the Grande detectors and the muonic component by the 
shielded array detectors of KASCADE. 

2. Test of hadronic interaction models 

The astrophysical interpretation of air shower data re- 
quires detailed knowledge of the properties of hadronic in- 
teractions at energies and kinematical ranges beyond the 
capabilities of present-day accelerator experiments. There- 
fore, air shower data are used to constrain hadronic inter- 
action models used in air shower simulations, such as the 
CORSIKA [10] code. Several hadronic interaction models 
have been systematically tested over the last decade. 

First quantitative tests [11112113] established QGSJET 98 
[TJ] as the most compatible code. Similar conclusions have 
been drawn for the successor code QGSJET 01 [15]. Re- 
cently, a new method has been developed to measure the 
attenuation of hadrons in air showers [16] . This method 
is very sensitive to inelastic hadronic cross sections. The 
investigations indicate, that the inelastic cross sections in 
QGSJET 01 are slightly too large (« 5% at 10 6 GeV). 
Predictions of QGSJET II [17|18|19] exhibit problems 



when compared to air shower data [20 , in particular, the 
predicted electron-hadron correlations are not compatible 
with the measurements. 

Predictions of SIBYLL 1.6 [21] were not compatible 
with air shower data, in particular there were strong in- 
consistencies for hadron-muon correlations. These findings 
stimulated the development of SIBYLL 2.1 [25]. This 
model proved to be very successful, the predictions of this 
code are fully compatible with KASCADE air shower data 
[23124115] , 

Investigations of the VENUS [25] model revealed some 
inconsistencies in hadron-electron correlations [13] . The 
predictions of NeXus 2 [26] were found to be incompati- 
ble with the KASCADE data, in particular, when hadron- 
electron correlations have been investigated [15] . Recently, 
predictions of the interaction model EPOS 1.61 [27|28|29j 
have been compared to KASCADE air shower data [3D] . 
The analysis indicates that EPOS 1.61 delivers not enough 
hadronic energy to the observation level and the energy per 
hadron seems to be too small. These findings stimulated the 
development of a new version EPOS 1.9 J3Tj. Correspond- 
ing investigations with this new version are under way and 
results are expected to be published soon. 

Analyses of the predictions of the DPMJET model yield 
significant problems in particular for hadron-muon corre- 
lations for the version DPMJET 2.5 [32], while the newer 
version DPMJET 2.55 is found to be compatible with air 
shower data [T5] . 

Presently, the most compatible predictions are obtained 
from the models QGSJET 01 and SIBYLL 2.1. 

3. Muons in air showers 

The lateral distribution of muons (E^ > 230 MeV) is in- 
vestigated using the shielded detectors of the original KAS- 
CADE array up to distances to the shower axis exceeding 
600 m [33] . A function suggested by Lagutin and Raikin is 
used to described the measured muon densities as a func- 
tion of the distance to the shower axis. This method allows 
to reconstruct the total number of muons in air showers 
with an accuracy better than 20%. 

A comparison of measured muon densities to predictions 
of air shower simulations indicates that these observable is 
sensitive to differences in hadronic interaction models [34 . 
The muon density has been measured as a function of the 
total number of electrons in the showers, i.e. as a function 
of primary energy. 

Muons with energies exceeding 800 MeV are investigated 
with the muon tracking detector. Lateral distributions have 
been measured, they can be described by a Lagutin-likc 
function [36J. The information of the reconstructed muon 
tracks is used to obtain radial and tangential angles with 
respect to the shower axis. These quantities, in turn, are ap- 
plied to infer the pseudorapidity distributions of secondary 
particles in hadronic interactions which decay into muons 
[37]. 
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Fig. 2. Effective depth of the muon production H£ in the N e — 
plane for showers with zenith angle smaller than 18° 1351 . Diagrams 
are overlayed for separate KASCADE and KASCADE-Grande anal- 
yses. 

The muon tracks are also used to measure the production 
height of muons in the atmosphere by means of triangu- 
lation [35] . The corresponding production depth of muons 
in the atmosphere is depicted in Fig. [2] in the N e — 
plane. The figure combines two separate analyses for the en- 
ergy ranges of KASCADE and KASCADE-Grande, respec- 
tively In the electron-muon number plane light particles are 
expected to be to the upper left of the main diagonal. In this 
region large muon production depths are reconstructed. 
This indicates light particles, penetrating deep into the at- 
mosphere. The measured air showers are classified in four 
intervals of the muon production depth. For these intervals 
the primary energy spectrum is reconstructed, based on the 
measured number of electrons and muons. This yields en- 
ergy spectra for four elemental groups. The spectra exhibit 
a fall-off for the light component at low energies (several 
10 15 eV), while the all-particle spectrum at high energies 
is dominated by heavy nuclei [35] . 

4. All-particle energy spectrum 

Different methods are applied to the data to reconstruct 
the all-particle energy spectrum. 

The density of charged particles at large distances from 
the shower axis is a good mass-independent estimator for 
the shower energy [38 . The particular optimum radial dis- 
tance depends on the layout and altitude of the detector ar- 
ray, for KASCADE-Grande a distance of 500 m yields best 
values and we use the charged-particle density S*(500) at 
500 m from the shower axis as energy estimator. The con- 
stant intensity cut method is applied to evaluate the atten- 
uation behavior of the S*(500) values. Thus, the all-particle 
energy spectrum is reconstructed. 

Another analysis takes advantage of the total number of 
charged particles measured for each air shower [3j5]. The 
shower size spectra measured in different zenith angle in- 
tervals are used to obtain the attenuation length for this 



observable. The constant intensity cut method is again ap- 
plied to correct for attenuation effects and the primary en- 
ergy spectrum is obtained. 

In a similar way the muon data have been analyzed |40] . 
The constant intensity cut method is used to correct for 
the attenuation of the measured total number of muons 
in different zenith angle intervals. Again, predictions of air 
shower simulations are used to establish an energy scale for 
the observed quantities and the all-particle energy spec- 
trum is derived from the observed number of muons. 

The last analysis discussed here, combines information 
from charged particles and muons [41] . Independent fits to 
the lateral distributions of charged particles and muons al- 
low to extract the number of charged particles and number 
of muons for each shower. These quantities are used to as- 
sign an energy to each shower on an event-by-event basis. 
Thus, the all-particle energy spectrum is inferred. 

The results obtained with the different methods are com- 
piled in Fig. [3] The reconstructed all-particle flux is shown 
as a function of the primary energy. The arrow indicates 
a typical systematic uncertainty of the energy scale of the 
order of 20%. Since the elemental composition is presently 
not clear in this energy region, for some methods two spec- 
tra are plotted, one for proton induced showers and a sec- 
ond one for iron induced showers. The different methods 
agree well for a heavy dominated mass composition around 
energies of 10 17 eV [43]. The spectra obtained in the en- 
ergy range between 10 16 and 10 18 eV smoothly extend the 
energy spectrum as measured by KASCADE to high ener- 
gies. The new results match as well the flux obtained by 
other experiments, shown in the figure as well. 

5. Composition 

First analyses have started to reveal the elemental com- 
position at the end of the galactic component [44]. They 
use the measured number of electrons and muons and in- 
dicate that these observables are indeed sensitive to the el- 
emental composition. At present, the measured data can 
be described by a fit of three elemental groups. A compar- 
ison of the measured electron-to-muon number ratio from 
KASCADE-Grande to similar data obtained with KAS- 
CADE yields consistent results. 

Ultimate goal is a deconvolution of the measured size 
spectra for the electromagnetic/charged and the muonic 
shower components, similar to the unfolding procedure ap- 
plied to KASCADE data [1] to obtain energy spectra for 
separate elemental groups at energies up to 10 18 eV. 

6. Conclusions and outlook 

The recent results from the KASCADE-Grande experi- 
ment indicate that we make substantial progress in mea- 
suring the energy spectrum and the elemental composi- 
tion of cosmic rays in the energy region between 10 16 and 
10 18 eV. The new data will significantly contribute to the 
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Fig. 3. All-particle energy spectrum of cosmic rays. Spectra obtained with different methods applied to KASCADE-Grande data are compared 
to results of other experiments I42I . 



understanding of the end of the galactic component and the 
transition to an extra-galactic component. The data will be 
crucial to distinguish between different astrophysical sce- 
narios. 
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